Influence of Mo on the Fe:Mo:C nano-catalyst thermodynamics for single- walled 

carbon nanotube growth 
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We explore the role of Mo in Fe:Mo nanocatalyst thermodynamics for low-temperature chemical 
vapor deposition growth of single walled carbon nanotubes (SWCNTs). By using the size-pressure 
approximation and ab initio modeling, we prove that for both Fe-rich (~ 80% Fe or more) and 
Mo-rich (~ 50% Mo or more) Fe:Mo clusters, the presence of carbon in the cluster causes nucle- 
ation of M02C. This enhances the activity of the particle since it releases Fe, which is initially 
bound in a stable Fe:Mo phase, so that it can catalyze SWCNT growth. Furthermore, the presence 
of small concentrations of Mo reduce the lower size limit of low-temperature steady-state growth 
from ^0.5SnmforpureFeparticlesto ~ 0.52nm. Our ab imtiothermodynamic modeling explains 
experimental results and establishes a new direction to search for better catalysts. 
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Critical factors for the efficient growth of single 
walled carbon nanotubes (SWCNTs) via catalytic chem- 
ical vapor deposition (CCVD) 1 '-! 3 - are the composi- 
tions of the interacting species (feedstock, catalyst, 
support^££££&i£), the preparation of the catalysts, and 
the synthesis conditions 1 ^^^^^' 1 -^. Efficient cat- 
alysts must have long active lifetimes (with respect to 
feedstock dissociation and nanotube growth), high se- 
lectivity and be less prone to contaminatio n 19 i 20 i 21 i 22 i 23 . 
Common factors that lead to reduction in catalytic activ- 
ity are deactivation (i.e. chemical poisoning or coating 
with carbon), thermal sintering (e.g. caused by highly 
exothermic reactions on the clusters surfac o 13 i 14 ' 24 with 
insufficient heat transfe r 25 i 26 ) and solid-state reactions 
(nucleation of inactive phases in the cluste r 22 i 23 ' 25 ). 

Metal alloy catalysts, such as Fe:Co, Co:Mo and 
Fe:Mo, improve the growth of CNTs4»iP-22i2L22i22 J 22 1 2I., 
because the presence of more than one metal species can 
significantly enhance the activity of a catalys t 20 ' 32 ' 33 ' 34 , 
and can prevent catalyst particle aggregatio n 20 ' 31 ' 32 ' 33 . 
In the case of Fe:Mo nanoparticles supported on AI2O3 
substrates, the enhanced catalyst activity has been 
shown to be larger than the linear combination of the in- 
dividual Fe/Al 2 3 and Mo/Al 2 3 activities' 2 -^. This 
is explained in terms of substantial inter-metallic inter- 
action between Mo, Fe and C 20 ' 35 ' 36 which is congruent 
with previously observed solid-state reactions between 
these elements. In fact, the addition of Mo in mechanical 
alloying of powder Fe and C mixtures 3 ^ promotes solid 
state reactions even at low Mo concentrations by forming 
ternary phases, such as the (Fe,Mo)23C6 type carbides 3 ^. 

The way in which carbon interacts with transition met- 
als depends on the metal species. Fe and Co belong to 
the "carbon dissolution-precipitation mechanism" group, 
where relatively large fractions of carbon dissolve into 



the cluster before stable carbides are formed, while Mo 
belongs to the "carbide formation-decomposition" group 
where carbide formation occurs rapidly at low carbon 
concentrations 3 ^. These mechanisms are governed by the 
interplay between solubility of C in the metal matrix (in- 
terstitial and/or substitutional defects) and the ease of 
metal carbide formation. Although relating C solubil- 
ity and catalytic ability of metal catalysts is not simple 
(factors like temperature, diffusion, kinetics and parti- 
cle size have to be considere d 22 ' 23 ), metals which form 
carbides can be efficient catalysts as long as they offer 
sufficient carbon solubility, even if the formation of the 
carbides interferes with the SWCNT growth, as is the 
case of Fo 22 ' 23 . Co can dissolve a relatively large fraction 
of carbon C and does not form stable ordered carbides, 
whereas Mo has limited C solubility and forms several 
stable carbides (even at low temperature). Similarly to 
Co, Fe can dissolve a relatively large fraction of C (at 
least in the 7 phase), but concomitantly it forms carbides 
at temperatures relevant to SWCNT growth (Fe 3 C)&22.. 
Nucleation of Fe 3 C and subsequent loss of C solubility 20 
can be used to identify the thermodynamic lower limit 
for the CVD growth of very thin SWCNTs from Fe 
nanoparticle o 22 ' 23 . Thus, for efficient CVD growth from 
Fe, Mo, Fe:Mo and Co:Mo clusters, the particles must be 
small (dispersed) and metallic (preserved from carbidiza- 
tion). In fact, active catalyst species in Co:Mo catalysis 
are also small metallic specie s 32 ' 33 with Mo added to limit 
particle aggregation by forming molybdate species which 
are later reduced to give small metallic Co:Mo nanopar- 
ticles. It has been shown that disruption of the Co-Mo 
interactions leads to a lower SWCNT growth efficiency^ -. 

The use of alloys instead of pure metals provides addi- 
tional degrees of freedom, such as fractions of the metal 
species, which has chemical and thermodynamic advan- 
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tages. The chemical advantages arise since the fractions 
of constituent species can be tailored to enhance catalytic 
performance that can be measured by the yield and qual- 
ity of nanotube o 32 ' 33 ' 34 . It has been found that low Fe:Mo 
ratios are favorable for growing SWCNTs (on AI2O3 sub- 
strates) since the presence, after activation, of the phase 
Fe2(Mo04)3 can lead to the formation of small metallic 
clusters^ 4 - (the best compromise between catalytic activ- 
ity and SWCNT selectivity was found to be Fe3.5M011.5j 
Refill). High fractions of Fe in Fe:Mo lead to formation 
of larger particles during the reduction step which are in- 
active for SWCNT growth, unless precautions are taken 
to avoid excessive sintering 2 ^ (no formation of molyb- 
date species has been reported when using CH4 feed- 
stock with Fe:Mo catalysts^ - and hence the chemical role 
of Mo in Fe:Mo is different from that in Co:Mo). The 
thermodynamic advantages are revealed when consider- 
ing the vapor- liquid-solid model (VLS), which is the most 
probable mechanism for CNT growt h 12 i 20 ' 41 . The metal- 
lic nanoparticles are very efficient catalysts when they 
are in the liquid or viscous states^ - , probably because 
one has considerable carbon bulk-diffusion in this phase 
(compared to surface or sub-surface diffusion). Gener- 
ally, unless stable intermetallic compounds form, alloy- 
ing metals reduce the melting point below those of the 
constituent s 35 ! 36 . This happens, for instance, with the 
addition of small fractions of Mo to Fe. Hence, to en- 
hance the yield and quality of nanotubes, one can tailor 
the composition of the catalyst particle to move its liq- 
uidus line below the synthesis temperature 2 ^. However, 
identifying the perfect alloy composition is non trivial. In 
fact, the presence of more than two metallic species al- 
lows for the possibility of different carbon pollution mech- 
anisms by allowing thermodynamic promotion of ternary 
carbides. So far, due to the very complex interplay be- 
tween competing phases at the nano-scale (even for pure 
Fe particle o 22 ' 23 ), the search for the proper catalyst com- 
position has been empirica l 20 ! 32 ' 33 ! 34 ' 40 . 

In the present manuscript we address the interaction 
between C and Mo:Fe nanoparticles. We use thermo- 
dynamic and quantum mechanical results to discuss the 
complexity of Fe- (~80% Fe or more) and Mo-rich (~50% 
Mo or more) Fe:Mo catalysts, previously addressed phe- 
nomenologically. The results are useful for the devel- 
opment of efficient catalysts for nanotube and graphene 
growth. 

Methods. Investigating the behavior of C in Fe:Mo 
nanoparticles requires an understanding of the interplay 
of the various phases of the Fe-Mo-C system at the nano- 
scale. Determining the thermodynamic stability of differ- 
ent phases in nanoparticles of different sizes with ah ini- 
tio calculations is complicated and computationally ex- 
pensive. In Re b 22 ' 23 we have developed a simple model, 
called the "size-pressure approximation" , which allows 
one to estimate the phase diagram at the nanoscale start- 
ing from bulk calculations under pressure. 



The size-pressure approximation. Surface curvature 
and superficial dangling bonds on nanoparticles are re- 
sponsible for internal stress fields which modify the 
atomic bond lengths. For spherical clusters, the phe- 
nomenon can be modeled with the Young-Laplace equa- 
tion AP — 2j/R where the parameter 7 (surface tension 
for liquid particles) can be calculated with ah initio meth- 
ods. As a first approximation, where all surface effects 
that are not included in the curvature are neglected, the 
study of phase diagrams for spherical particles can be 
mapped onto the study of phase diagrams for bulk sys- 
tems under the same pressure that is produced by the 
curvature. It is important to mention that 7 is not a 
real surface tension but an ah initio fitting parameter de- 
scribing size-induced stress in nanoparticles. In our case, 
it can be assumed that 7 is independent of the fraction 
of Mo (since the best Fe:Mo catalysts are Fe-rich) and 
C (since the amount of C in the catalytically active par- 
ticle state is limited) 10 ' 20 . Figure 2 of Ref32 shows the 
implementation of the "size-pressure approximation" for 
Fe nanoparticles. The idea is simple. Two interpolations 
are involved: "pressure versus distortion" and "distor- 
tion versus curvature" are coupled to obtain a relation 
"pressure versus curvature" (where curvature is 1/i? 22 ). 
In more details, the left hand side of Figure 2 of Ref. 22 
shows the average distortion of the bond length inside 
the cluster Ad nn = d° n — d nn for a variety of spherical 
bec particles as a function of the inverse radius (\/R). 
The right hand side shows the compression of the bond 
length as a function of the hydrostatic pressure for the 
bulk system. Combining the two sets of data yields: 

P- R= 2.46 GPa-nm (7 = 1.23J/m 2 ). (1) 

Equation ([1} is used to deduce the Fe-Mo-C phase di- 
agram of nanoparticles of radius R from ah initio cal- 
culations of the bulk material under pressure P. The 
parameter 7 compares well with the experimental values 
of the surface tension of bulk Fe at the melting point ~ 
1.85 J/m 2 — and ^1.90^. For a detailed explanation of 
the "size-pressure approximation" see Reb2£. 

The assumption of 7 being independent of the Mo con- 
centration is justified at zero (and low) temperature be- 
cause molybdenum tends to segregate inside the particle, 
so that it can not affect the chemistry and the bonding 
states at the surface (see the section "surface energies" 
near the end of the manuscript). At high temperature, 
Mo would eventually populate the surface and modify 
7. To a first approximation, by linearly interpolating 
between the experimental values of the surface tension 
for Fe (1.85 J/m 2 ) and Mo (2.08 J/m 2 ) (i.e. Vegard's 
law), and assuming that the parameter 7 follows the same 
trend, we would obtain a ~ 3% increase for 7 at the 
optimal composition Fe4Mo of the catalyst (see below). 
Thus, all the estimated radii and diameters deduced from 
equation [T] might be overestimated of few percent at high 
temperature (energies and pressures are not affected by 
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7)- Considering the formidable computational needs re- 
quired to obtain a concentration-dependent parameteri- 
zation of 7, we believe that our calculated value is accu- 
rate enough to describe the general trends of size-induced 
thermodynamics and activity. 

Quantum mechanical calculations. Simulations are 
performed with VASP 44 , using projector augmented 
waves (PAW)^ 5 - and exchange-correlation functionals as 
parameterized by Perdew, Burke and Ernzerhof (PBE) 46 
for the generalized gradient approximation (GGA). Sim- 
ulations are carried out with spin polarization, at zero 
temperature, and without zero-point motion. All struc- 
tures are fully relaxed. Numerical convergence to within 
about 2 me V/ atom is ensured by enforcing a high en- 
ergy cut-off (500 eV) and dense k-meshes. The hydro- 
static pressure estimated from the pressure-size model is 
implemented as Pulay stress 4 ^. Ternary phase diagrams 
are calculated using bcc-Mo, bcc-Fe and SWCNTs as ref- 
erences (pure-Fe phase is taken to be bec because our 
simulations are aimed at the low temperature regime of 
catalytic growth). The reference SWCNTs have the same 
diameter of the particle to minimize the curvature-strain 
energy. In fact, CVD experiments of SWCNT growth 
from small (^0.6-2.1 nm) particles indicate that the di- 
ameter of the nanotube is similar to the diameter of the 
catalyst particle from which it grows. In some exper- 
iments where the growth mechanism is thought to be 
root-growth, the ratio of the catalyst particle diameter 
to SWCNT diameter is ~1.0, whereas in experiments in- 
volving pre-made floating catalyst particles this ratio is 
~1.6~. Formation energies are calculated with respect to 
decomposition into the nearby stable elements or phases, 
depending the position in the ternary phase diagram as 
described later. 

Competing phases. Pure elements (bcc-Fe, bcc-Mo and 
SWCNT-C) are included in the calculation in their most 
stable low temperature configurations. Other phases are 
included if they are stable in the temperature range used 
in CVD growth of SWCNTs or if they have been reported 
experimentally during or after the growt h 35 i 36 i 37 . Thus, 
we include the binaries Mo 2 C, Fe 2 Mo and Fe 3 G2 5 .' 3 £ 
In addition, since our Fe-rich Mo:Fe experiments were 
performed with compositions close to Fe4Mo2°-, we in- 
clude a random phase Fe4Mo generated with the special 
quasi-random structure formalism (SQS). Bulk ternary 
carbides, which have been widely investigated due to 
their importance in alloys and steel, can be considered 
as derivatives of binary structures with extra C atoms 
in the interstices of the basic metal alloy structures. 
Three possible ternary phases have been reported for 
bulk Fe-Mo-C 4 ^ and they are referred as t\ (MqC), t 2 
(M 3 C) and t 3 (M23C6) (M is the metal species). For 
simplicity, we follow the same nomenclature. t\ is the 
well-known MqC phase, which has been observed ex- 
perimentally as Fe4Mo 2 C and Fe 3 Mo 3 C structures (77 
carbides) 50 ' 51 . Both of these structures are fco 4 ^ but 



have different lattice spacings 5 ^. Our calculations show 
that the most stable variant n is Fe4Mo2C, and we de- 
note it as T\ henceforth. r 2 is the Fe2MoC phase, which 
has an orthorhombic symmetry distinct from that of 
FeaC 5 ^.. We consider FeziMozCgS 5 . as the third t 3 fee 
phase. We use the Cr 23 Cg as the prototype structure 5 ^ 
where Fe and Mo substitute for Cr. Although M 23 Cg 
type phases do not appear in the stable C-Fe or C- 
Mo systems, they have been reported in ternary C-Fe- 
Mo systems and also appear as transitional products 
in solid state reactions^ 9 -. Time-temperature precipita- 
tion diagrams of low-C steels have identified r 2 -Af 3 C, 
t 3 -M 23 C6 and t\-MqG as low-temperature, metastablc 
and stable carbides, respectively^. Furthermore, r 2 - 
M 3 C carbides precipitate quickly due to carbon-diffusion 
controlled reaction while r 3 -M 23 Cg carbides precipitate 
due to substitutional-diffusion controlled reactions. The 
latter phenomenon, requiring high temperature, longer 
times and producing metastable phases 5 ^ is not expected 
to enhance the catalytic deactivation of the nanoparti- 
cle. In summary, as long as the presence of carbon does 
not lead to excessive formation of M 3 C (Fe 3 C and r 2 - 
Fe 2 MoC), the catalyst should remain catalytically active 
for SWCNT growth. 

Results. A structure at a given composition is consid- 
ered stable (at zero temperature and without zero-point 
motion) if it has the lowest formation energy of all struc- 
tures at this composition and, if on the ternary phase 
diagram, it lies below the convex hull of tie planes con- 
necting all the other stable structure a 57 ' 58 ' 59 ' 60 . Phases 
lying above the convex hull and with small positive for- 
mation energies may be favored at higher temperatures 
due to configurational and vibrational entropy contribu- 
tions. 

We generate the convex hulls with the qconvex 
package®- 1 . By projecting the 3D facets onto the 2D 
plane we obtain the graphs shown in Figure [TJ where 
the panels (a),(b),(c) and (d) represent the phase di- 
agrams at zero temperature of nanoparticles of radii 
R ~ oo,1.23,0.62,0.41 nm, calculated at P = 0,2,4, 
and 6 GPa, respectively. Stable and unstable phases are 
shown as black squares and red dots, respectively. The 
solid green lines, connecting the stable phases are the 
projected edges of the 3D convex hull "facets" . The num- 
bers 1...8 in panels (c) and (d) indicate the boundary 
phases' crossings of our test cases, and are discussed be- 
low. Table U shows the formation energies of the various 
competing phases calculated with respect decomposition 
into the appropriate stable components (reported in the 
"refs." lines). 

In each panel of Figure [TJ the two dotted lines con- 
necting Fe 4 Mo and FeMo to SWCNT denote the intro- 
duction of carbon into the system. Fe4Mo has been re- 
ported to be an effective catalyst composition^ while 
FeMo represents a hypothetical Fe:Mo particle with a 
Mo content larger than 33%. This Mo-rich nanoparti- 
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FIG. 1: (color online). Ternary phase diagram for Fe-Mo-C 
nanoparticles of R ~ oo, 1.23, 0.62, 0.41 nm. Notations are 
explained in the text. 

cle gives low yields of SWCNTs at high temperature^. 
In thermodynamic terms, the growth of SWCNTs has 
to be considered as nucleation of ordered C, hence the 
balance between all the competing phases along the dot- 
ted lines controls the activity of the particle. The size 
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TABLE I: Formation energies (meV/atom) for binary and 
ternary phases for nanocatalysts of different sizes: FeaC, 
Mo 2 C, Fe 2 Mo, n, r 2 and Tz, calculated with respect to the 
reference species shown in the table. 



of the particles have been calculated with the interpo- 
lation shown in Figure 2 of Refj 22 ' 23 . We investigate a 
very small particle with N atoms ~ 40, which may able 
to support SWCNT growth^ 3 -, to explore size-induced 
stabilization trends. As seen from the zero temperature 
phase diagrams in Figure [T] and the energies in Table [H 
M02C and Fe2Mo are stable for bulk materials as well 
as nanoparticles, while FeaC and T3 are stable only for 
small nanoparticles of radius R ~ 0.41 nm (E/[Fe3C]< 
for R < 0.58 n m 22 i 23 ). t\ and T2 are always unstable. 

Fe^Mo particles. The most obvious advantage that a 
Fe4Mo particle has over a pure Fe particle is that the 
[Fe4/ 5 Mo 1 / 5 ]i_ a: -C a ; line does not intercept any carbide 
(FesC, T3,T2). This implies that, at least at low temper- 
atures, there is a surplus of unbounded metal (probably 
even at high temperatures since the line is far from all 
of the competing stable phases). This is illustrated in 
Figure [3J which shows the fractional evolution of species 
as one progresses along the [Fe4/ 5 Mo 1 / 5 ]i_ a ;-C 2 ; line in 
Figure [T] 

Large Fe^Mo particles (R > 0.62 nm). When a large 
Fe4Mo particle {R > 0.62 nm in Figured] (a)) is exposed 
to carbon feedstock, the M02C phase nucleates by disso- 
ciating Fe2Mo. This is seen for concentrations between 
< x c < 0.09. The vertical dashed line "1" in the figure 
indicates the boundary crossing of the [Fe4/ 5 Mo 1 / 5 ]i_ 2; - 
C x pollution path with the Fe^Mo2C edge, as shown in 
Figure QJc). Nucleation of M02C has two consequences: 
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FIG. 2: (color online). Panel a). Fraction of species for cata- 
lyst composition [Fe 4 / 5 Mo 1/ / 5 ]i_ 2 :-Ca; and nanocatalyst of R > 
0.62 rim. The dashed vertical line, labeled as "1", represents 
[Fe4/ 5 Mo 1 /5]i_ I -C I crossing the boundary phase Fe«Mo2C, 
as shown in Figure [j7c). Panel b). Fraction of species for 
catalyst composition [Fe4 / / 5 Mo 1 / 5 ]i_ :E -C I and nanocatalyst 
of R < 0.41 nm. The dashed vertical lines, labeled as "4" 
and "5", represent [Fe4 / / 5 Mo 1 /5]i_ :c -Ca; crossing the boundary 
phases Fe<->Mo2C and T3«Mo2C, as shown in Figure QJd). 



it consumes carbon atoms that are added to the particle 
by bonding them to molybdenum, and it releases free-Fe 
which is beneficial for the catalytic activity of the clus- 
ter (the Fe fraction increases with C concentration in the 
range < x c < 0.09). In addition, while forming M02C, 
carbon does not produce SWCNTs. In fact, we need to 
saturate the particle in C and enter the region Fe-Mo2C- 
SWCNT of the phase diagram in Figure QJc) before free 
carbon is available. Hence, the fraction of C is non-zero 
only after x c ~ 0.09 in Figure Ufa). For x c < 0.09 the 
free Fe is expected to be on the particle surface since 
the fraction of free-Fe is larger than that of M02C and, 
as presented later, surface energy calculations show that 
free-Fe resides at the surface of the particle whereas the 
Mo (M02C) parties found near particle core. To con- 



clude, steady state growth of SWCNTs is possible from 
large Fe4Mo particles since free, catalytically active Fe is 
present, even for x c > 0.09. 

Small Fe^Mo particles (R < 0.41 nm). Similarly to 
the large Fe4Mo particle, initial exposure of the small 
Fe4Mo cluster (i? < 0.41 nm) to carbon feedstock pro- 
motes M02C nucleation by dissociating Fe2Mo and re- 
leasing free Fe. However, for x c > 0.09 the scenario is 
different. In fact, while the fraction of free Fe increases 
between < x c < 0.09 (left of line "4" in Figured] (b)), 
it reduces to zero between 0.09 < x c < 0.23 (points "4" 
and "5" of Figure Bid)). Concomitantly, the carbide T3 
nucleates by bonding carbon atoms that are added to 
the particle and hence eliminating the amount of free Fe 
that was released during the growth of M02C. Similarly 
to Fe-C particle a 22 ' 23 with the carbide FesC, the nucle- 
ation of T3 is clearly detrimental for SWCNT growth. 
In fact, the lack of simultaneous presence of free-Fe and 
excess C causes the particle to be catalytically inactive. 
However, the different formation energies of T3 and Fe3 
lead to different size thresholds between the two regimes 
(growth/no-growth). In Refs i 22 ' 23 we estimated that the 
minimum radius needed for a low-temperature Fe par- 
ticle to be active is R^ in ~ 0.58 nm. For the Fe4Mo 
case, by interpolating the energies of T3 reported in Ta- 
ble HI by determining P at which E[ts(P)] = 0, and by 
using the size/pressure approximation relation {1}, we 
obtain R^ Mo ~ 0.52 nm. The reduced size advantage 

^rnin Mo < -^min agrees with experimental results that 
show that nanotubes, grown from Fe-rich Fe:Mo catalysts 
are thinner than those grown from pure Fe particles 2 ^. As 
pointed out befor e 22 ' 23 , the presence of i?^^ M ° defines 
three possible thermodynamic scenarios: steady-state-, 
limited- and no-growth for R > R%? n Mo , R ~ R^ n Mo 
and R < JJ^-* M °, respectively. It is worth mentioning 
that as long as the concentration of Mo in a Fe-rich Fe:Mo 
clusters is chosen to be slightly larger than that of T3 we 
obtain ^„ M ° < R^ in . 

FeMo test particles. To address the poor growth capa- 
bility of Mo-rich Fe:Mo catalysts (low yield of SWCNTs 
at high temperature^ 2 -) we repeated the above analysis 
with an equiconcentration FeMo particle. For simplicity, 
we choose to avoid dealing with the Mos.iFeoy.g phase 
reported at low temperatur e 35 ' 36 because, as our results 
show, its presence does not affect the pollution mecha- 
nism discussion. 

Large FeMo particles (R > 0.62 nm). Contrary 
to Fe4Mo particles, large FeMo clusters contain ample 
amounts of molybdenum capable of nucleating M02C. Ir- 
respective of whether Mos.iFeoy.g is included in the dis- 
cussion, both free-Fe and an excess of C are present for 
x c > 0.20 (line "3" in Figure (3Ja) and phase boundary 
"3" in Figure QJc)). Hence, the growth SWCNTs is pos- 
sible for x c > 0.20. However, since the fraction of Fe 
considerably smaller than that of M02C, the activity of 
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FIG. 3: (color online). Panel a). Fraction of species for 
catalyst composition [FeirtMoi/ali-x-Qr and nanocatalyst 
of R > 0.62 nm. The dashed vertical lines, labeled as "2" 
and "3", represent [Fex^Moi^i-i-Ca: crossing the bound- 
ary phases Fe2Mo<->Mo2C and Fe<->Mo2C, as shown in Fig- 
ure OIc). Panel b). Fraction of species for catalyst com- 
position [FeirtMoi^i-a-Cx and nanocatalyst of R < 0.41 
nm. The dashed vertical lines, labeled as "6", "7" and "8", 
represent [Fei^Moi^i-x-Ca, crossing the boundary phases 
Fe2Mo<->Mo2C, Fe^Mo2C, and T3^->Mo2C, as shown in Fig- 
ure Did). 



the whole FeMo particle is drastically smaller than that 
of the Fe4Mo catalyst. Thus, although it is possible to 
growth SWCNTs from large FeMo clusters, the expected 
yield is low and the synthesis temperature needs to be 
high (to overcome the reduced fraction of catalytically 
active free-Fe), as experimentally reported 62 . Thermo- 
dynamically, an excess of C and a amount of free-Fe for 
x c > 0.20 guarantees the existence of steady-state growth 
of SWCNTs (albeit slow and inefficient). 

Small FeMo particles (R < 0.41 nm). Small FeMo par- 



ticles are similar to small Fe4Mo clusters. Nucleation of 
T3 and the absence of free-Fe and excess C means that 
the particle is catalytically inactive. The minimum clus- 
ter radius, R^^° is the same as i?^ M ° because this 
quantity is determined by the stabilization of the same 
T3 phase. Thus, there are three possible thermodynamic 
scenarios that are similar to the previous case, i.e., low- 
yield steady-state-, low-yield limited- and no-growth for 
R > Bg£°, R ~ R^ n 10 and R < R™f°, respectively. 
The analysis can be extended to different fractions of Fe 
and Mo and the results are summarized in Table HT1 



catalyst 


Fe-rich Fei-^Mo^ 


Mo-rich Fei-^Mo^ 


x c range 


x* 3 < x < 1/3 


x > 1/3 


(test) 


Fe 4 Mo 


FeMo 


p ^ TjFeMo 


steady-state growth 


low-yield steady-state growth 


R ~ R FeMo 


limited growth 


low-yield limited growth 


p . pFeMo 


no growth 


no growth 



TABLE II: Scenarios of thermodynamically allowed SWCNT 
growth modes for different sizes and compositions of Fe:Mo 
catalysts. x* 3 is the concentration of Mo in T3 removed of all 
C. R^t'° = 0.52 < = 0.58 nm from Refs.-^^. 



Surface energies. The chemical species that are 
found on the surface of a Fe:Mo cluster with coexisting 
Fe, and Fe2Mo phases is the species that has the lowest 
surface energy per unit area, 7. This quantity can be 
obtained from the relation 7 = [E n — nEs] /2A, where 
E n is the total energy of a n-layer slab, Eb is the total 
energy of a single bulk layer, and A is the surface of the 
unit cell (the factor two accounts for the creation of two 
surfaces) 64 . For bcc Fe, we construct slabs with the low- 
est surface energy (1 1 0) termination 64 . Surface termina- 
tion for the ordered Fe2Mo phase has not been reported in 
experiments or in theoretical works. By using the pack- 
age AFLOW, which performs high-throughput simultane- 
ous optimization of planar density (high) and number of 
broken bonds (low) 64 i 66 , we find that the closed packed 
(0 4) plane of Fe2Mo has the lowest energy. We obtain 
7 (! 10 ) = 2.44 J/m 2 and j^ Mo = 3.12 J/m 2 . The rela- 
tion 7i? e < 7Fe 2 Mo indicates that Mo will not be at the 
cluster surface. Thus, in Fe-rich Fe:Mo nanocatalysts, Fe 
covers as much surface area as possible, and the afore- 
mentioned nucleation of M02C caused by C pollution is 
advantageous for the activity of the particle, by releasing 
free-Fe on the surface. 

Conclusions. In this paper we discuss the role of Mo 
in the thermodynamic properties of Fe:Mo nanocataly- 
sists by using the size-pressure approximation and ah ini- 
tio modeling. We show that for both Fe-rich (~ 80% Fe 
or more) and Mo-rich (~ 50% Mo or more) Fe:Mo clus- 
ters, the presence of carbon causes nucleation of M02C 
which enhances the activity of the particle by releasing 
free-Fe. With respect to pure Fe-catalysts, the addition 
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of Mo (up to small concentrations) decreases the size of 
the smallest catalyst needed for low-temperature CVD 
steady-state growth to R^^f° ~ 0.52 nm. 
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